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We introduce dimensional perturbation techniques to 
Bose-Einstein condensation of inhomogeneous alkali gases 
(BEC). The perturbation parameter is 5 = where k de- 
pends on the effective dimensionality of the condensate and on 
the angular momentum quantum number. We derive a sim- 
ple approximation that is more accurate and flexible than the 
N oo Thomas-Fermi ground state approximation (TFA) of 
the Gross-Pitaevskii equation. The approximation presented 
here is well-suited for calculating properties of states in three 
dimensions and in low effective dimensionality, such as vortex 
states in a highly anisotropic trap. 



I. INTRODUCTION 

The most commonly used approach to describe a di- 
lute gas of atoms in a BEC at T = is mean-field the- 
ory, which takes the form of the time-independent Gross- 
Pitaevskii equation (GPE): 

= M^(r), (1) 

where the three-dimensional coupling constant is U3 = 
(47r?i^a)/m, a is the s-wave scattering length, Vtrap is the 
external trapping potential, fi is the chemical potential, 
and N is the number of condensate atoms. The com- 
plex order parameter, '0(r), is referred to as "the wave 
function of the condensate." 

The N —^ 00 Thomas-Fermi approximation (TFA) has 
been proven to be a highly successful analytical approx- 
imation of the GPE @J|l. The strength of the cx3 
TFA is its simplicity: neglecting the kinetic energy re- 
sults in a simple approximation of the ground-state con- 
densate density that is effective in analyzing properties 
of large-N condensates. For condensates with a mod- 
erate number of atoms and condensates with attractive 
interactions, the TFA breaks down. The kinetic energy 
is important in each case, especially in the latter, where 
the kinetic energy is necessary to prevent collapse. The 
effects of attractive interactions have been studied using 



approximation techniques such as variational trial wave 
functions |^-|^. Other approximations that have been 
employed to extend the TF regime of validity include the 
?i — > TFA 1^, a method that uses two-point Pade ap- 
proximants between the weakly- and strongly-interacting 
limits of the ground-state p|, and a variational method 
for anisotropic condensates 1 1] . 

With the study of BEC in highly anisotropic traps as 
a motivation, we use a perturbation formalism that per- 
mits the effective dimensionality, D, to vary. Because 
it readily allows one to approximate quantities in any 
dimension, such a formalism is ideally suited for con- 
densates in the anisotropic traps used in many laborato- 
ries where the condensate can be effectively one-, two- 
or three-dimensional. The perturbation parameter is 
S = where k depends on the effective dimension- 

ality of the condensate and on the angular momentum 
quantum number. The <5 ^ limit becomes an exactly 
soluble problem, the solution of which is used by the vari- 
ous dimensional-scaling methods as the starting point for 
the solution of the full three-dimensional problem . 
The (5 — *■ approximation to the condensate density, 
which retains part of the kinetic energy, is quite accu- 
rate for both a large and moderate number of atoms in 
the BEC ground state, and the dimensional scaling for- 
malism, which treats the dimensionality as a parameter, 
is advantageous when studying condensates of low effec- 
tive dimensionality due to extreme trap anisotropy. The 
centrifugal term in the (5 — > density also makes it a 
good physical starting point for treating vortex states. 



II. iV ^ 00 THOMAS-FERMI APPROXIMATION 
TO THE GROUND-STATE 

In the case of positive scattering length, the repulsive 
interaction causes the density to become flat, and the 
kinetic energy of the condensate becomes negligible in 
the N ^ 00 limit. This limit of the GPE results in the 
highly successful classical approximation for the density 
of the condensate ground-state known as the Thomas- 
Fermi approximation {N ^ 00 TFA). The N ^ 00 TFA 
for the ground state in a three-dimensional isotropic trap 

is 01 



PTF{r) = 



1 



NU: 



1 



TF 



muj r ) 



(2) 



*e-mail: brett_mckinney@mailaps.org 



e-mail: dwatson@ou.edu 



for /xtf > 



1 9 9 



and ptf — elsewhere. Eq. (0) pro- 



vides an excellent description of the condensate ground- 
state density in the bulk interior. This approximation 
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breaks down near the surface of the gas where the den- 
sity is not flat; the wave function must vanish smoothly, 
making the kinetic energy appreciable in the boundary 
layer. The chemical potential is obtained from continu- 
ity and normalization of Eq. (^. In oscillator units, one 
finds 



(3) 



where R = {IBNa/ahoY^^ is the Thomas-Fermi classical 
cutoff radius in oscillator units of length, aho = \/nJrriuj. 

Boundary layer theory techniques have been employed 
to obtain corrections to the N ^ oo TFA at the conden- 
sate surface where the gradient of the density is no longer 
small p"^-p5| . The leading order correction to the ground 
state chemical potential due to the boundary layer at the 
surface is of order R~^ln{R). 



III. EFFECTIVE DIMENSIONALITY 



et. al For A ^ 1, one may assume the motion in 

the radial-direction in the x-y plane is described by a 
gaussian of small width, x(7'_l), and, following the same 
procedure, one obtains a ID equation analogous to Eq. 
(Q), where Ni would represent the number of atoms in 
the ID condensate per unit area in the xy-plane. 



IV. GPE IN VARIABLE DIMENSIONALITY 

We begin by explicitly generalizing the nonlinear 
Schrodinger equation (NLSE), Eq. (^, to D-dimensions 
where r becomes the radius of a D-dimensional sphere 
with D-1 remaining angles. The Laplacian is general 



ized to D-dimensions (Bohn, Esry and Greene [g8 19 
treat the Laplacian in a similar fashion, in which they 
use hyperspherical coordinates to define a mean conden- 
sate radius.), and the potential terms retain their three- 
dimensional form; the coupling constant is generalized in 
the final scaling. We obtain the Schrodinger equation: 



We use a perturbation formalism where the effective 
dimensionality, D, of the condensate is allowed to vary. 
The effective-dimensionality of the condensate depends 
on the relative size of the confinement length in each of 
the three spatial dimensions. Most experimentally real- 
ized traps are axially symmetric with some having a high 
degree of anisotropy |16|Jl7|] . In the case of axial symme- 



try, the trapping potential takes the form, Vtrap{f) 
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-\- A^z^), where A — 



LOzjiox^ is a measure of the degree of anisotropy. The sys- 
tem reduces to a three-dimensional isotropic condensate 
for A = 1. In the small- (large-) A limits, the system 
reduces to an effective one- (two-) dimensional isotropic 
condensate. It is conceivable that an isotropic hamilto- 
nian in a fractional-dimensional space could be used to 
describe experimental condensates for intervening values 
of A, but we will focus our attention on integer dimen- 
sions. 

As an illustration, consider A 3> 1 where the motion of 
the atoms in the z-direction becomes frozen and their mo- 
tion is described by a gaussian of small width. To deter- 
mine the 2D effective coupling constant, we assume the 
wave function in Eq. (|^) is separable: ■(/;(r) = V'2(''±)x(-2^)j 
where x[^) is assumed to be a gaussian, and operating 
with J dzx' 1 one finds a new effective 2D GPE: 
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which has the same form as Eq. ([l|), but r is the 2D radius, 
= = M - ^Wz/2, N2 = N J dz\x\^- Requiring 

that ■02 and x be normalized to unity, / has units 

of 1 /length; thus, we interpret N2 as the number of atoms 
in the 2D condensate per unit length along the z-axis. 
In our subsequent scalings, we will adopt a notation for 
the number of atoms that is similar to that of Jackson 
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+7V;73|V'(r)|'|^(r)=A^^(r), (5) 

where Lj~,_^ is a generalized angular momentum operator 
depending on D-1 angles with eigenvalues —1{D + 1 — 2) 
|20| ; the angular momentum quantum number, I, is non- 
negative. Substituting these eigenvalues and introduc- 
ing the radial Jacobian factor in a transformation of the 
wave function, 0(r) — A-'^^^'>^'^ip{r), to eliminate the first 
derivative terms, we find 
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Finally, we make two sets of scalings to arrive at the 
NLSE in dimensionally scaled oscillator units. The first 
scaling is a purely dimensional scaling: r = K^f, Q = k^uj, 
jl = K^yit, and = K^0, where k = D + 21. The final 
scaling is to scaled oscillator units (denoted by bars) : r = 



o-hoT, A — hoj^i, and ip 
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Combining these two scalings, we arrive at 



tp, where a^o — ^JhJmCj. 
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where everything is now in dimensionally scaled oscilla- 
tor units and 8 = 1/k. For the effective dimensions of 
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primary interest in this study, the dimensionally scaled 
coupling constants are, for 3D, §3 = g^/K^^"^, where 
93 = AirN^a/aho and N3 is the number of condensate 
atoms; and for 2D, g2 — §2/^^, where 52 = 47rA^2a and 
N2 represents the number of atoms in the 2D condensate 
per unit length along the z-axis. The definition of N2 
makes (72 dimensionless. For general D, 
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making Eq. valid for describing condensates in any 
effective dimension. In the next section, we describe a 
simple and accurate zeroth-order approximation to Eq. 
(0). 



where n{D) = 27r^/^ /r{D/2^ The 5 ^ limit can be 
thought of as a large-D or large-1 limit. 

Eq. (||) is valid where the density is non-negative. In 
addition to the N ^ 00 TF-like classical cutoff radius 
near the surface, Rmax, the centrifugal term requires that 
another cutoff be defined, Ro, slightly removed from the 
origin, to satisfy the requirement that the density be non- 
negative. In terms of the chemical potential, the cutoff 
radii in scaled oscillator units are defined as 



and 
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In regular oscillator units (a^o). 
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V. ZEROTH-ORDER DENSITY 

It has been pointed out by Schuck and Vifias that 
the true TF limit (7i ^ as originally applied to the 
case of Fermi statistics |2l|]) is not equivalent to iV cxd, 
and they show that the ?i ^ TF limit for bosons does 
not neglect the kinetic energy for the ground state. The 
N ^ 00 TEA to the ground state is too harsh on the 
kinetic energy for a moderate number of atoms. A less 
harsh and nearly as simple approximation is the zeroth- 
order ((5 0) approximation of Eq. (0). Unlike the 
ground-state N ^ 00 TFA, which neglects the entire ki- 
netic energy, our zeroth-order approximation of the gen- 
eralized GPE neglects the derivative part of the kinetic 
energy but retains a centrifugal term. For vortex states, 
one understands this term as being a centrifugal barrier 
due to quantized circulation, which pushes atoms away 
from the axis of rotation. This centrifugal barrier arises 
from the condensate phase: jV^p {tp = ^/pe''^, where p 
is the condensate density and S is the spatially dependent 
condensate phase; then the condensate velocity is given 
by V = ^VS"). For the ground state, this centrifugal 
term in the zeroth-order density has an alternate, quan- 
tum mechanical interpretation, which helps explain its 
good agreement with numerical calculations. We discuss 
this interpretation in Section [v| . 

The 6 —^ limit of the angular-dimensional perturba- 
tion parameter in Eq. results in the following zeroth- 
order density in scaled oscillator units: 



go 8r^ 2 



for Ro{p) < r < RmaxifJ-) and p — elsewhere, 
normalization condition becomes 
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Notice for the ground state in the strongly interacting 
regime that p ^ 1 and the cutoff radii for the ground 
state become N ^ 00 TF-like: i?o ~ and p « R^^^./2\ 
the strongly interacting limit, or, equivalently, the N 
00 limit of our zeroth-order approximation collapses to 
the N ^ 00 TFA, as expected. (For finite N, as will be 
shown later, our zeroth-order approximation gives better 
agreement with the numerical solution of the GPE than 
the N ^ 00 TFA.) Using the integration limits defined 
in Eq. (pi), along with the condensate density defined 
in Eq. (|^y7the normalization condition (Eq. [l0| ) gives an 
equation for the zeroth-order chemical potential that is 
easily solved in any dimension. (See Section VII where 
this procedure is illustrated for two dimensions.) 

Once the chemical potential is calculated, it is then 
used in Eq. (^ to complete the description of the zeroth- 
order wave function. One can then calculate the energy 
from 



E/N = J d 
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^As the GPE is nonlinear, one cannot treat excited-/ states 
(vortices) for D > 2 radially symmetric traps in the usual 
manner of separating the wave function into radial and an- 
gular parts. Presently, however, vortices in 3D radially sym- 
metric traps are not realized. Vortex states in a 2D isotropic 
trap do not pose a problem to theory because the spheri- 
cal harmonic wave function acts as a phase factor. In ID, 
the spherical harmonic wave function is a constant and, since 
there are no angles, one can think of / = and I = 1 as even 
and odd parity states. 
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or in the zeroth-order approximation and scaled oscillator 
units, 



E/N^n{D) J r 

« Emu/N + Eho/N + E„,t/N. 
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VI. GROUND STATE IN THREE DIMENSIONS 

The numerical effect of the centrifugal-like term in the 
zeroth-order approximation on the ground state of a sta- 
tionary condensate is clear from Fig. |l|, where we compare 
the numerical solution of the GPE chemical potential 
with our zeroth-order approximation and the N oo 
TFA for up to 10, 000 atoms in a spherical trap. 

Our zeroth-order approximation is more accurate than 
the iV — > oo TFA for all N, most notably for a moderate 
number of atoms. The accuracy of the zeroth-order ap- 
proximation is comparable to boundary layer corrections: 
the zeroth-order approximation is slightly more accurate 
for small coupling constant, while boundary layer theory 
is slightly more accurate for larger coupling constant, but 
the difference between all three approximations becomes 
small for very large coupling constant. 

The correct physical interpretation of this centrifugal- 
like term, as originally noted by Chatterjce [Q, is that it 
is the component of the kinetic energy needed to satisfy 
the minimum uncertainty principle. The zeroth-order 
density includes a centrifugal term from the kinetic en- 
ergy, which pushes the wave function away from the ori- 
gin in the ground state as if there were a non-zero quan- 
tum of angular momentum; however, the 1/r^ contribu- 
tion to the ground state density of a nonrotating cloud 
clearly is not due to any rotational motion of the cloud. 
This effect, which becomes less pronounced as N in- 
creases, is demonstrated in Figs. ^ and |[ which show the 
numerically calculated GPE ground-state non-Jacobian 
weighted wave function (^) along with our zeroth-order 
approximation and the N oo TFA. 

The centrifugal-like term in the lowest order of dimen- 
sional perturbation theory {6 —^ 0) can be understood 
as arising from the requirement that the system's uncer- 
tainty product be a minimum |^^. Another way to see 
how a centrifugal term may arise in the ground state - 
this time within the N oo TFA - is by applying the 
Langer modification of WKB theory to the N ^ oo TF 
density. For vortices, one may not neglect the entire ki- 
netic energy in the iV — > oo limit. A slightly more general 
N oo TF density than Eq. that includes vortices 
is 
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For a spherical trap, A^ = l{l + 1), which reduces to the 
usual ground state — > oo TF density for / = 0, but 
using the Langer modification, where the correct asymp- 
totic phase of the WKB wave function is obtained by the 
replacement l{l -I- 1) ^ (/ -I- 1/2)^ in the centrifugal po- 
tential, a centrifugal barrier remains in the ground-state: 
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The dependence of our perturbation parameter on the 
angular momentum quantum number suggests that the 
zeroth-order density will be a good physical starting 
point for vortex states, which we explore in the next sec- 
tion for D = 2. The remaining centrifugal term in our 
zeroth-order approximation is a lowest-order correction 
to the kinetic energy, which, for a moderate number of 
atoms, greatly improves the ground state approximation 
{I — 0) over the — > oo TFA. For a very large number of 
atoms, the contribution from the kinetic energy becomes 
very small, as can be seen by comparing the wave func- 
tions in Figs. I and |, for A^ = 10"^ and A^ 10^ ^'^Rb 
atoms, respectively. As the number of atoms increases, 
our unphysical core becomes smaller than the healing 
length, eventually vanishing: our zeroth-order and the 
N oo TF wave functions become indistinguishable 
from the numerical solution for large N. 



VII. LOWER DIMENSION 

The i5 ^ density (Eq. |9|) is well suited for describ- 
ing condensates in the presence of a vortex, where the 
centrifugal term models the vortex core (see Fig. ^). 
In this section, we present explicit expressions for the 
D — 2 ground-state and vortex states. In the angular- 
dimensional scaling of the GPE, one has considerable 
freedom in the choice of the scaling parameter, 6 — 1/k. 
In the previous section, we used k — D + 21 or, for the 
ground state, k = D. The choice k = D + 21 — 2 ex- 
actly reduces our expressions below for the chemical po- 
tential and energy to a two-dimensional N —f oo TFA 
that includes the leading contribution to the kinetic en- 
ergy due to fluid motion of the condensate . Slightly 
improved agreement of the zeroth-order energy with the 
numerical solution of the 2D GPE for a wide range of 
mean-field coupling constant can be obtained by choos- 
ing K — D + 21 — 1, which changes the numerator in the 
centrifugal term of Eq. (|^) to 1 — Zeroth order predic- 
tions show a small amount of variability with the choice 
of K, but the results of higher order perturbation theory 
should not depend on the particular choice of k. 

Using Eqs.(|jl^Jll|) for D = 2, one finds that the scaled 
chemical potential satisfies 



2^=2V'^ 4+16^" 7 



(17) 
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Recalling the earlier conversion relations leading to Eq. 
(0), the chemical potential, in regular oscillator units, 
satisfies 



92 / 2 ^ ^2 / M 



(18) 



Solving Eq. (^7|) for the scaled zeroth-order chemical po- 
tential, /2, and using the resulting wave function, Eq. (||) 
and Eq. (|l|) give a simple analytical approximation for 
the 2D energy, 
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where we have already converted to regular oscillator 
units. Eq. ( |T^ ) and Eq. ( p^ for the chemical potential 
and energy per atom, respectively, are analagous to the 
results of the N ^ oo TEA for vortices given in Rcf. , 
which includes a kinetic energy term associated with the 
fluid motion that is encoded in the wave function's phase. 
This similarity is due to the zeroth-order 5—^0 limit be- 
ing a large angular momentum limit or, in the language 
of hydrodynamics, a large quantum of circulation limit. 
The zeroth-order approximation for D = 2 results in a 
shifted TF-like energy spectrum, whose ground state ap- 
proximation is, just as for D — 3, more accurate than the 
N ^ oo TEA for any coupling constant, most noticeably 
for smaller coupling. For the energy of a single charge 
vortex located at the center of the trap, the above expres- 
sions are more accurate than the unregulated N oo 
TEA in Ref. |Q for a moderately sized coupling con- 
stant, and slightly less accurate for very large coupling. 
Eor our zeroth-order approximation and N oo TEA, 
respectively, the relative errors in the first vortex state 
energy are 0.56% and -0.88% for 52 = 1, 000; and 0.015% 
and -0.012% for c/2 = 100, 000. 



for the chemical potential. This improvement is caused 
by the centrifugal-like term, which brings in the kinetic 
energy needed to satisfy the minimum uncertainty prin- 
ciple and which adds a needed outward push to the wave 
function. The accuracy of the zeroth-order approxima- 
tion is comparable to the lowest order correction due to 
the boundary layer at the condensate surface. Improved 
accuracy for the ground state is most noticeable for a 
moderate number of atoms, the case in which the ki- 
netic energy is most significant. For a sufficiently large 
number of atoms with positive scattering length, the ki- 
netic energy becomes small for the ground state chemical 
potential, and the three approximations converge to the 
numerical solution of the GPE. 

The core near the origin and the presence of the angu- 
lar momentum quantum number, I, in the scaling param- 
eter, 6 — 1/k, make the zeroth-order (6 0) density an 
especially good starting point for studying properties of 
vortices. The ground-state N ^ 00 TEA is unable to ac- 
commodate such states, but it can be extended to include 
vortices by introducing the gradient of the phase from the 
Laplacian [p3|-p6[. We expect higher order, finite-i5 cor- 
rections to further refine the shape of our zeroth-order 
density for the ground state and vortex states. 

We have shown that the dimensional scaling formalism 
is conducive to analysis of condensates of any dimension. 
We outlined how simple yet accurate approximations can 
be achieved for any effective D, and we demonstrated the 
improved numerical results for D = 3. In addition to 
3D BEG, the dimensional scaling formalism provides a 
useful analytical tool in the study of BEG in lower effec- 
tive dimensionality. In future work, we plan to test the 
approximation presented here on other observables and 
states of BEG in D-dimensions, and we are extending the 
methods presented herein to D-dimensional cylindrical 
coordinates, where the anisotropy parameter is included 
explicitly for treatment of axially symmetric traps with 
arbitrary anisotropy. 



VIII. CONCLUSIONS 
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FIG. 1. Chemical potential in oscillator units vs. number 
of condensate atoms for a *^_R6 condensate in a 3D isotropic 
trap, where a = 100 bohr and 1/ — 200 Hz. The zeroth-order 
{6 — > 0) approximation of dimensional perturbation the- 
ory presented here (dashed) is in better agreement with the 
numerical solution of the GPE (solid) than the N —> 00 
Thomas- Fermi approximation (dash-dot). 
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FIG. 3. Same as Fig. |^ but with 10® atoms, corresponding 
to gs 87, 204. As the number of atoms increases, the un- 
physical core in our zeroth-order wave function shrinks. Near 
the origin, the N ^ 00 TF and numerically calculated wave 
functions overlap, while the N ^ 00 TF and our zeroth-order 
wave function overlap in the boundary region. For sufficiently 
large N, the three wave functions become indistinguishable. 



numericai 

- N ->» 
- - S->0 




radius (osc. units) 



FIG. 2. Ground state wave functions {tp - non-Jacobian 
weighted) for a Rb condensate of 10,000 atoms in a 3D 
isotropic trap, where a = 100 bohr and v = 200 Hz. These pa- 
rameters correspond to gs ~ 872.04. Plotted are the numer- 
ical solution of the GPE (solid), the N 00 Thomas- Fermi 
approximation (dash-dot) and our zeroth-order (5 — > 0) ap- 
proximation (dashed). Our zeroth-order approximation con- 
tains an unphysical core near the origin, but the added kinetic 
energy, which causes the core to appear, is also responsible for 
the increased accuracy seen in Fig. 0. 
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FIG. 4. Comparison of the condensate wave function (-0 - 
non-Jacobian weighted) with anl — 1 vortex in a 2D isotropic 
trap with g2 — 10, 000. The solid line is the numerical solution 
of the GPE and the dashed line is our zeroth-order {S 0) 
wave function, whose centrifugal term models the vortex core. 
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